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Abstract 

Exclusive p + p~ production in two-photon collisions between a quasi-real photon, 
7, and a mid-virtuality photon, 7*, is studied with data collected at LEP at centre- 
of-mass energies 183 GeV < ^/s < 209 GeV with a total integrated luminosity of 
684.8 pb _1 . The cross section of the 77* — > p + p~ process is determined as a function 
of the photon virtuality, Q 2 , and the two-photon centre-of-mass energy, W 77 , in the 
kinematic region: 0.2 GeV 2 < Q 2 < 0.85 GeV 2 and 1.1 GeV < W jy < 3 GeV. These 
results, together with previous L3 measurements of p°p° and p + p _ production, allow 
a study of the 77* — > pp process over the Q 2 -region 0.2 GeV 2 < Q 2 < 30 GeV 2 . 
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1 Introduction 



The L3 Collaboration has recently measured the exclusive production of p°p° [1,2] and p + p [3] 
pairs in the two-photon fusion process: 

e + e~ — > e + e~77* — > e + e~pp , (1) 

where the beam electrons 1 -* radiate virtual photons which interact and produce a hadronic final 
state. One of the photons, 7, is quasi-real, characterised by a small value of its squared four 
momentum, P 2 = m 2 ~ 0, whereas the other one, 7*, has a significant virtuality, Q 2 = — P 7 » 2 = 
— m 7 » 2 — m 2 . Our measurements cover the two-photon centre-of-mass energy 

1.1 GeV < W, n < 3 GeV. (2) 

The two measurements [1,3] done at large virtualities, 1.2 GeV 2 < Q 2 < 30 GeV 2 , provide 
a testing ground for a recently-developed QCD-based model [4]. This model describes well the 
Q 2 -dependence of the p°p° production at large momentum transfer [5]. The measured cross 
sections for p°p° and p + p~ production were found to have a similar dependence on and 
to be of similar magnitude. However, the p + p~ cross section is systematically higher than 
the p°p° one. This is in contrast with the suppression and different W iy dependence of p + p~ 
production [6] with respect to p°p° [7] observed in data with Q 2 pa and < 2 GeV. We 
note that despite the wide range of theoretical models [8,9], p-pair production at Q 2 pa is still 
not well understood. Therefore the experimental study of the (Revolution of p-pair production 
is important to understand vector meson pair-production in two-photon interactions. 

Previously, we performed a measurement of p°p° production [2] for intermediate virtualities: 

0.2 GeV 2 < Q 2 < 0.85 GeV 2 . (3) 

In this Letter, we complement that study with the first measurement of the process: 

e + e~ — > e + e~77* — > e + e~p + p~ (4) 

in the kinematic region (2) and (3). These data allow to follow the (Revolution of the pp- 
production over two orders of magnitude in this variable. 

The analysis techniques employed in this study are similar to those of our previous mea- 
surements [2,3]. The data used, corresponding to an integrated luminosity of 684.8 pb _1 , are 
the same as in Reference 2 and were collected by the L3 detector [10] at LEP at centre-of-mass 
energies 183 GeV < ^/s < 209 GeV. Scattered beam electrons which have radiated photons 
with virtualities in the range (3) can be "tagged" by the Very Small Angle Tagger (VSAT) [11]. 
The VSAT is an electromagnetic calorimeter, constructed with BGO crystals, with a geomet- 
rical acceptance covering the polar angle range 5 mrad < 8 < 10 mrad, for azimuthal angles in 
the ranges —1.25 rad <0< 1.25 rad and 7r— 1.25 rad <0<7r+1.25 rad. When the electron with 
the largest scattering angle is detected in the VSAT, the virtuality of the photon it radiated is, 
within 1% precision, equal to the transverse momentum squared, p 2 , of the final state hadron 
system: 

Q 2 = 2E b E s {\ - cos0 s ) w E b E s 9 2 Pa p 2 , (5) 

where E b is the beam energy, and E s and 6 S are the energy and the scattering angle of the 
tagged electron. Therefore the VSAT is not used to directly measure Q 2 , but rather to select 
exclusive final states by correlating the direction of the transverse momentum vector of the 
tagged electron with the detected hadron system. 

^ Throughout this Letter, the term "electron" denotes both electron and positron. 
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2 Event Selection 



The reaction e + e 



e + e p + p contributing to the process 




(6) 



is identified by one and only one scattered electron, e tag , detected in the VSAT, two charged 
pions measured in the tracking chamber, and energy clusters from the two-photon decays of 
the 7r°'s, deposited in the BGO electromagnetic calorimeter. These events are collected by two 
independent track-triggers [12]. The trigger efficiency, as determined from the data itself, is 



Single-tagged events are selected by requiring just one electromagnetic cluster with energy 
greater then 50% of the beam energy reconstructed in the VSAT. The event candidates must 
have exactly two tracks with zero total charge. The tracks must come from the interaction 
vertex, have transverse momentum greater than 100 MeV and an energy loss in the tracking 
chamber compatible with the pion hypothesis. The selected events should contain a 7r°7r° pair, 
therefore we consider event candidates that have four or five photons, identified as isolated 
clusters in the electromagnetic calorimeter. Photons having energies greater than 60 MeV are 
paired to reconstruct neutral pions, which are required to be in the mass window 100 MeV < 
■^(77) < 170 MeV, as shown in Figure la. The mass of a tt° candidate is constrained to the 
nominal value by a 1-C kinematic fit. If more than one 7r°7r° combination exists, the one with 
the smallest x 2 sum °f the fits is taken. To make the selection robust against instrumental noise 
and backgrounds and to reduce the sensitivity to the Monte Carlo simulation of fake photons, 
we retain events with one additional photon, not used in the 7r°7r° pair, if the photon energy is 
less than 300 MeV and does not exceed 10% of the energy of the 7r°7r° pair. 

The transverse momentum squared, p 2 , of the four-pion system is used to measure the Q 2 
of the event and is required to be in the range 0.2 GeV 2 — 0.85 GeV 2 . For selection of an 
exclusive final state, the acoplanarity angle, aco , calculated from the difference between the 
azimuthal angle of the tagged electron, (ptag, shown in Figure lb, and the azimuthal angle of 
the four-pion system, is required to be less than 0.4 rad, as shown in Figure lc. The data 
contain a contribution from rj production, as visible in the the 7r + 7r~7r° mass spectrum, shown 
in Figure Id. This background is removed by requiring M(7r + 7r~7r°) > 0.65 GeV. 

After all cuts, 414 events are retained. Their four-pion mass spectrum is shown in Figure 2a. 
The region 1.1 GeV < W 77 < 3 GeV is populated by 387 events, which are used for the 
cross section determination. A strong signal from production is observed in the M(7r ± 7r°) 
spectrum, shown in Figure 2b. The clustering of entries at the crossing of the mass bands in 
the correlation plot of the masses of the 7i ± n combinations, shown in Figure 2c, gives evidence 
for a signal from p + p~ intermediate states. No structure is observed in the correlation plot 
of the masses of the ir + ir~ and -k°-k° combinations, shown in Figure 2d. We also inspected 
the two- and three-pion mass distributions, shown in Figure 3, for production of higher-mass 
resonances. The only statistically-significant signal is from the af (1320) state in the tt ± -k tt 
mass spectrum, as seen in Figure 3f. 



(60 ± 3)%. 
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3 Data Analysis 



3.1 Monte Carlo Modelling 



To estimate the number of p + p events in the selected four-pion data sample, we consider 
non-interfering contributions from the processes: 



About 40 million Monte Carlo events of the processes (7) are generated with the EGPC [13] 
program, which uses the luminosity function from Reference 14. Particle production and decay 
is uniform in phase-space. The generated events are passed through the full L3 detector simu- 
lation using the GEANT [15] and GHEISHA [16] programs and processed in the same way as 
the data, reproducing the detector behaviour as monitored in the different data-taking periods. 

For acceptance calculations, Monte Carlo events are assigned a Q 2 -dependent weight, eval- 
uated using the GVDM form-factor [17] for both interacting photons. The detection efficiencies 
of the process (4) are listed in Tables 1 and 2 for bins in Q 2 and W J7 . The efficiencies for the 
four-pion final states of all the processes (7) are of similar magnitude. 

3.2 Background Estimation 

The contribution to the selected events from e + e~ annihilation and from the process e + e~ — > 
e + e _ r + T - is negligible. Random coincidences with off-momentum beam electrons, which give 
signals in the VSAT, are a source of background. The flux of these particles is dominantly on 
the outer side of the LEP ring. Therefore, this background would cause an excess in the number 
of events having a tag on the outer side of the accelerator ring, N out , with respect to the inner 
side, N in . In the selected data, the ratio N out /N in = 1.04 ± 0.10 is close to unity, indicating 
that this background is small. This conclusion is corroborated by the good agreement observed 
between the (f> tag distribution of the selected data and Monte Carlo event samples, shown in 
Figure lb. 

Two sources of background remain. The first is partially-reconstructed events from two- 
photon interactions with higher particle multiplicities, when tracks or photons escape detection. 
The second is signal events with one or more photons substituted by photon candidates due to 
noise. To estimate the accepted background we use background-like event samples extracted 
from the experimental data. The first background is modelled with selected 7r ± 7r ± 7r 7r events, 
in which at least two charged particles have not been detected and by 7r + 7r~7r°7r 7r° events in 
which one n° is excluded from consideration. An event-mixing technique is employed in order 
to reproduce events from the second background: one or two photons forming a tt° are excluded 
from a selected event and replaced by photons from another data event. The <p aco distributions 
of the background-like data samples, passing the selection, are combined with the distribution 
of selected 7r + 7r~7r°7r° Monte Carlo events so as to reproduce the (f> aco distribution observed in 
the data, as shown in Figures lc. The estimated background levels are listed in Tables 1 and 2. 
As data samples are used in the background estimation, they contain also a fraction of events 
with fake tags and thus take into account the effect of this background. 



77 
77 
77 
77 




(7) 



4(1320)tt t ; 



7r + 7r 7i 7i , non— resonant. 
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3.3 Fit Method 



In order to determine the differential p + p~ production rate, a maximum likelihood fit of the 
data to a sum of Monte Carlo samples of the processes (7) is performed in intervals of Q 2 and 
W 77 using a box method [1-3,18]. The inputs to the fit are the six two-pion masses in an event, 
namely the four combinations 7r ± 7r and the two combinations 7r + 7r~ and 7r°7r°. They provide a 
complete description of a four-pion event in our model of isotropic production and phase space 
decay. 

The analysis procedure is optimised for deriving the p + p~ contribution and only the p + p~ 
content and the sum of the rest of the contributing processes, denoted as "other 4-ir" , are 
considered for cross section measurements. To check the quality of the fit, the two- and three- 
pion mass distributions of the data are compared in Figure 3 with those of a mixture of Monte 
Carlo event samples from the processes (7), in proportions determined by the fit. The observed 
experimental distributions are reasonably well described by the Monte Carlo model. 

4 Results 

The cross sections of the process e + e~ — > e + e~p + p~ in bins of Q 2 and W 77 , Aa ee , are listed 
in Tables 1 and 2. The statistical uncertainties, also listed in Tables 1 and 2, are those of the 
fit. The differential cross section, da ee /dQ 2 , derived from Aa ee , is listed in Table 1. When 
evaluating the differential cross section, a correction based on the Q 2 -dependence of the p + p~ 
Monte Carlo sample is applied, so as to assign the cross section value to the centre of the 
corresponding Q 2 -bin [19]. 

To evaluate the cross section, <7 77 , of the process 77* — > p + p~ , the integral of the transverse 
photon luminosity function, L TT , is computed for each Q 2 and W 7J bin using the program 
GALUGA [20], which performs 0(a 4 ) QED calculations. The same procedure was used in 
our previous studies [1-3]. The cross section a 77 is derived from the measured cross section 
using the relation <r 77 = A<j cc /Ltt- Thus, <r 77 represents an effective cross section containing 
contributions from both transverse and longitudinal photon polarisations. The cross section of 
the process 77* — > p + p~ is listed in Table 1 as a function of Q 2 and in Table 2 as a function of 
W 77 . The sum of the cross sections of the other contributing processes is also given in Tables 1 
and 2. 

Several sources of systematic uncertainty are considered. The contribution of the selection 
procedure is in the range 12% — 18%; Monte Carlo statistics in the range 1.3% — 2.1%; the 
fit procedure in the range 11% — 20%. Half of the changes of the acceptance when no form 
factor re-weighting of the Monte Carlo events is performed is considered as model uncertainty. 
It is in the range 0.5% — 5%. The background correction procedure introduces systematic 
uncertainties in the range 2% — 6%. All contributions are added in quadrature to obtain the 
systematic uncertainties, quoted in Tables 1 and 2. Finally, a normalization uncertainty of 5% 
accounts for the uncertainty of the trigger efficiency determination. 

5 Study of 77* — ► pp Process 

Combining the present results with the L3 data on pp production from References 1-3, we 
compare the p + p~ to the p°p° channels and their evolution as a function of Q 2 . The cross 
section of the process 77* — > pp is plotted in Figure 4 as a function of W 77 . For W 77 < 
2.1 GeV and 0.2 GeV 2 < Q 2 < 0.85 GeV 2 there is a clear enhancement of p°p° production 
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relative to p + p~ . This is similar to what was observed at Q 2 ~ [6, 7], but in contrast 
with the high Q 2 -region, where both cross sections have similar dependence on W yi and the 
p + p~ is systematically higher than the p°p°- These differences are clearly seen in the ratio 
R = £ Aa ee (p + p~ ) / E Ao- ee (p°p°) where the sum is for the region 1.1 GeV < W 77 < 2.1 GeV. 
In the domain 0.20 GeV 2 < Q 2 < 0.85 GeV 2 , we measure R = 0.62 ± 0.10 (stat.) ± 0.09 (syst.), 
a value that can only be explained by the presence of an isospin 1 = 2 intermediate state or 
by a mixture of different states [8,9]. The value of this ratio for 1.2 GeV 2 < Q 2 < 8.5 GeV 2 
is R = 1.81 ± 0.47 (stat.) ± 0.22 (syst.) [3], close to the factor 2, expected for an isospin 1 = 
state. Such variation suggests different p-pair production mechanisms at low and high Q 2 . 

The differential cross section da cc /dQ 2 of the reaction e + e~ — > e + e~pp is shown in Figure 5a. 
The L3 measurements span a Q 2 -region of two orders of magnitude, over which the differential 
cross sections show a monotonic fall of more than four orders of magnitude. The pp data are 
fitted to a form [21] expected from QCD-based calculations [22]: 

dff cc 1 (H) 

dQ 2 ~ Q n (Q 2 + < W jy > 2 ) 2 ' 1 ] 

where n is a constant and < > is the average 1V 77 value, 1.9 GeV for this measurement. 
Although this formula is expected to be valid only for Q 2 ^> W 77 , we find it provides a 
good parametrisation of the Q 2 -evolution of the pp data. A fit to the p + p~ data finds an 
exponent n = 2.3 ± 0.2 with x 2 /d.o.f. = 1.4/7. A value n — 2.9 ± 0.1 was found for p°p° with 
X 2 /d.o.f. = 6.9/10 [2]. Only the statistical uncertainties are considered in the fits. The results 
of the fits are shown in Figure 5a. The fits indicate a cross-over of the differential cross sections 
of p + p~ and p°p° production in the vicinity of Q 2 ps 1 GeV 2 . 

The measured cross section of the process 77* — > pp as a function of Q 2 is shown in 
Figure 5b. The change of the relative magnitude of p + p~ and p°p° production is clearly visible 
when comparing the low- and the high-Q 2 regions. A parametrisation, based on the GVDM 
model [17]: 

1 4- O 2 /Am 2 fl 99 

a^W^Q 2 ) = a^W iy )F(Q 2 ), and F(Q 2 ) = ^ r v (1 ± |^ + ± , (9) 

with r p = 0.65, r u = 0.08, = 0.05 and m = 1.4 GeV reproduces well the Q 2 -dependence 
of the p°p° data, as shown in Reference 2 and indicated by the line in Figure 5b. The fit 
finds a cross section of 13.6 ± 0.7 nb for the region 1.1 GeV < W yi < 3 GeV at Q 2 = 0. 
The (Revolution of p + p _ data cannot be satisfactorily described by this form. In addition, as 
shown in Figure 5b, the p°p° data cannot be described by the much steeper Q 2 -fall of a p-pole 
parametrisation [2]; the same is true for the p + p~ cross section since it is decreasing with Q 2 
less steeply than the p°p° one. 



6 Conclusions 

We have performed the first measurement of exclusive p + p~ production in mid-virtuality two- 
photon collisions. These results complement previous L3 measurements of exclusive p + p~ and 
p°p° production and allow to follow the evolution of pp cross sections over a Q 2 -region of two 
orders of magnitude. 

A QCD-based form, derived for the description of the differential cross-section da ee /dQ 2 
of the process e + e~ — > e + e _ pp at high Q 2 , is found to provide a good parametrisation of the 
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(Revolution of the pp data in the entire interval 0.2 GeV 2 < Q 2 < 30 GeV 2 , over which the 
differential cross sections show a monotonic decrease of more than four orders of magnitude, 
for 1.1 GeV < W 77 < 3 GeV. 

The Q 2 -dependence of the cross section of the process 77* — > p°p° is well reproduced by a 
parametrisation based on the GVDM model over the entire Q 2 -region. On the other hand, the 
p + p~ data cannot be satisfactorily described by such a parametrisation. A p-pole description 
of the Q 2 -dependence for both p°p° and p + p~ data is excluded. 

The relative magnitude of p + p~ and p°p° production changes in the vicinity of Q 2 ~ 1 GeV 2 , 
suggesting different p-pair production mechanisms at low and high Q 2 . 
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Q 2 range 




Bg 


Aa ee [pb] 


daJdQ 2 [pb/GeV 2 ] 


<t 77 [nb] 


<r 77 [nb] 


[GeV 2 ] 


[%] 


[%] 


p+p- 


p+p- 


p + p- 


other An 


0.20 - 0.28 


0.8 


14 


7.4 ±2.4 ±1.9 


92 ± 29 ± 23 


5.7 ± 1.8 ± 1.4 


10.9 ± 2.2 ± 1.5 


0.28 - 0.40 


1.2 


14 


5.7 ±1.8 ±1.3 


47 ±15 ±10 


4.3 ± 1.4 ± 1.0 


12.2 ± 1.8 ± 1.4 


0.40 - 0.55 


1.1 


15 


5.6 ±1.6 ±1.1 


37 ±11 ± 7.3 


4.9 ± 1.4 ± 1.0 


13.3 ± 2.0 ± 1.8 


0.55 - 0.85 


0.7 


18 


7.7 ±2.5 ±2.0 


25 ± 8.2 ± 6.5 


5.3 ± 1.7 ± 1.4 


12.1 ± 2.2 ± 1.9 



Table 1: Detection efficiencies, e, background fractions, Bg, and cross sections of the reactions e + e~ — > e + e~p + p~, 77* — > p + p~ and 
of the sum of the rest of the contributing processes, "other 4n" , as a function of Q 2 for 1.1 GeV < W 77 < 3 GeV. The values of the 
differential cross sections are corrected to the centre of each bin. The first uncertainties are statistical, the second systematic. An 
overall normalization uncertainty of 5% for the trigger is not included. 



W / 77 -range 




Bg 


Aa ee [pb] 


<r 77 [nb] 


<t 77 [nb] 


[GeV] 


[%} 


[%} 


p+p- 


p+p- 


other An 


1.10 - 1.40 


0.6 


25 


4.9 ± 1.8 ± 1.3 


3.9 ±1.5 ±1.1 


9.0 ± 2.4 ± 1.7 


1.40 - 1.65 


0.9 


18 


6.7 ± 1.6 ± 1.3 


7.6 ± 1.9 ±1.5 


14.8 ±2.7 ±2.5 


1.65 - 1.85 


1.1 


15 


5.1 ± 1.5 ±0.9 


8.4 ± 2.4 ± 1.6 


15.8 ±3.1 ±2.3 


1.85 - 2.10 


1.1 


13 


3.9 ± 1.4 ±0.8 


5.9 ±2.0 ±1.2 


18.3 ±3.0 ±2.7 


2.10 - 2.40 


1.2 


10 


2.2 ± 1.0 ±0.5 


3.2 ± 1.4 ±0.8 


11.5 ±2.1 ± 1.8 


2.40 - 3.00 


1.2 


10 


2.2 ± 1.0 ±0.5 


1.9 ±0.9 ±0.5 


8.5 ±1.5 ±1.5 



Table 2: Detection efficiencies, e, background fractions, Bg, and cross sections of the reactions 
e + e~ — > e + e~p + p~, 77* — > p + p~ and of the sum of the rest of the contributing processes, 
other 47r, as a function of W 77 for 0.2 GeV 2 < Q 2 < 0.85 GeV 2 . The first uncertainties are 
statistical, the second systematic. An overall normalization uncertainty of 5% for the trigger is 
not included. 
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Figure 1: Distributions for 7r + 7r~7r°7r° candidates, a) Two-photon invariant mass for the selected 
7r°'s (two entries per event); b) azimuthal angle, (fitag, of the tagged electron for tags in the inner 
side of the LEP ring (in) and, folded over it, for tags in the outer side of the LEP ring (out); c) 
acoplanarity angle, aco , between the tagged electron and the 7r + 7r~7r°7r° system and d) mass 
of the 7r + 7T~7r° system (two entries per event). The data are compared to the four-pion Monte 
Carlo. The estimated background is indicated by the hatched histograms. The arrows indicate 
the selection cuts. 
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Figure 2: Mass distributions for the selected events: a) the four-pion system, W 77 ; b) the 
combinations (four entries per event); c) correlation between the 7r~7T° and 7i + n° pairs (two 
entries per event) and d) correlation between the 7r + 7r~ and 7r°7r° pairs. The two-dimensional 
distributions have a bin width of 55 x 55 MeV 2 , the size of the boxes is proportional to the 
number of entries and both plots have the same vertical scale. 
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Figure 3: a),c),e),g) Mass distributions of the combinations (four entries per event) in 

four (^-intervals and distributions for the entire kinematic region 1.1 GeV < < 3 GeV and 
0.2 GeV 2 < Q 2 < 0.85 GeV 2 of b) The sum of the n+ir' and 7r°7r° mass spectra (two entries 
per event), d) The neutral three-pion combinations (two entries per event), f) The charged 
three-pion combinations (two entries per event), h) The sum of the 7r + 7r~7r° and tt^tt ?! mass 
spectra (four entries per event). The points represent the data, the hatched areas show the 
p + p~ component and the open areas show the sum of the other contributing processes. The 
fraction of the different components are determined by the fit and the total normalisation is to 
the number of the events. 
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Figure 4: Cross section of the process 77* — > pp as a function of W 77 in three Q 2 regions. The 
results from this measurement, full points in a), are compared to previous L3 measurements 
of the pp production [1-3]. The bars show the statistical uncertainties. Some points from the 
previous measurements are displaced horizontally for better readability. 
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Figure 5: The pp production cross section as a function of Q 2 , for 1.1 GeV < W 77 < 3 GeV: 
a) differential cross section of the process e + e~ — > e + e~pp and b) cross section of the process 
77* — > pp. The results from this measurement, full points in the region Q 2 < 1 GeV 2 , are 
presented together with previous L3 measurements of the pp production [1-3]. The bars indi- 
cate the statistical uncertainties. Some points from the previous measurements are displaced 
horizontally for better readability. The lines in a) represent the results of fits using the QCD- 
inspired form of equation (9). The lines in b) represent the results of a fit to the p°p° data 
based on the GVDM model [17] and of a fit based on a p-pole parametrisation. 
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